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1 Introduction 



The modern treatment of the pion-nucleon system as an effective field theory of the 
standard model was pioneered by Gasser, Sainio, Svarc [|l|] and Krause who formulated 
the "relativistic" version of baryon chiral perturbation theory. It was then shown by 
Jenkins and Manohar that the methods of heavy quark effective theory allow for 
a systematic low-energy expansion of baryonic Green functions in complete analogy to 
the meson sector. The latter approach is usually called heavy baryon chiral perturbation 
theory. 

Applications of this effective field theory beyond the tree level require the knowledge 
of the divergences generated by one-loop graphs. For the pion-nucleon interaction in the 
heavy mass expansion, the full list of one-loop divergences to 0{p^) has been worked out 
by Ecker 0. This analysis was then extended to the three- flavour case by Miiller and 
MeiBner 0. In these papers, the bosonic loop and the mixed loop (boson and fermion 
lines in the loop) were treated separately. This required a cumbersome investigation of 
the singular behaviour of products of propagators, because the mixed loop does not have 
the form of a determinant, like the purely bosonic or fermionic loops. 

To overcome these difficulties, we have recently developed a method where bosons 
and fermions are treated on an equal footing. Employing the notions of supermatrices, 
superdeterminants and supertraces H, 0, we have constructed a super-heat-kernel rep- 
resentation for the one-loop functional of a boson-fermion system. In this approach, the 
determination of the one-loop divergences is reduced to simple matrix manipulations, in 
complete analogy to the familiar heat-kernel expansion technique for bosonic or fermionic 
loops. 

The present paper is organized as follows: in Sect. |^ I briefly review the super- heat- 
kernel method. In contrast to the Euclidean space formulation used in [0] , the presentation 
in this work refers to Minkowski space throughout. In Sect. ^ the super-heat- kernel for- 
malism is applied to a rather general class of scalar-heavy fermion interactions (including, 
of course, heavy baryon chiral perturbation theory). The one-loop divergences to second 
order in the fermion fields are given explicitly. These results are then specialized to the 
two-flavour version of heavy baryon chiral perturbation theory in Sect. ^. My conclusions, 
together with an outlook to possible extensions of the present work, are summarized in 
Sect. ^. Several momentum-space integrals are collected in the Appendix. 

2 Super-Heat-Kernel 

Let us consider a general action 

S[ip,ilj,i^]= J d'^x C{^,ij,i^) (2.1) 

for ub real scalar fields ipi and np spin 1/2 fields ipa- Anticipating the later use of dimen- 
sional regularization, I am starting in d-dimensional Minkowski space. To construct the 
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generating functional Z of Green functions, these fields are coupled to external sources 
ji (i = 1, . . . , ns), Pa, Pa ia = l,...,nF), 

Z[j,p,p] = e^^[J>'^] = / [dyDdijdiJi e^(5[^'V',V^]+j^^+^;p+PV') ^ (2.2) 



where W[j,p,p] is the generating functional of connected Green functions. I have used 
the notation 

f(fi + ^p + -p^:= / d'^X {jiipi + IpaPa + Pai^a) ■ (2.3) 



The normalization of the functional integral is determined by the condition Z[0, 0, 0] = 1. 
The solutions of the classical equations of motion 

dLfi dip a dlpa 

are denoted by ipd, ipd- They are uniquely determined functionals of the external sources. 
With fluctuation fields ^, r] defined by 

V'a = V'cl.a + Va , (2.5) 



the integrand in (|2.2| ) is expanded in terms of ^, rj, rj. The resulting loop expansion of the 
generating functional 

W = Wl=o + Wl=i + ... 
starts with the classical action in the presence of external sources: 

Wl=o = S[(pd, i>ch ^ci] + j'^Vci + -ipciP + y#ci • (2.6) 
The one-loop term Wl=i is given by a Gaussian functional integral 

e'^^=' = f [d^d7]dT]] e^^"^[^^"^^"^ci;€,r;,»j] ^ ^2.7) 



where 



is quadratic in the fluctuation variables. Employing the notation introduced in (|2.3|) , S^"^^ 
takes the general form 
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= l(e^Ae + ?ji?r7-r7^5V + e^rr7-r/^r^e + ^re-e^rV) , (2.9) 

where A,B,r,r are operators in the respective spaces; A = and B are bosonic 
differential operators, whereas F and F are fermionic (Grassmann) operators. They all 
depend on the classical solutions ipch i^ci- 
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The standard procedure for the evaluation of (|2.7| ) is to integrate first over the fermion 
fields r], r] to yield the bosonic functional integral 



This leads to the familiar result 



L=l 



lndet(v4 - VB-'V + V'^B-^^V^) - Indet Aq - z(lndet B - Indet 5, 



- Tr In — 
2 



i Tr In 



B 

B'n 



Tr ln(l - A-^TB-^T + A-^T^B-^^T^) 



i ^ , A „ B ^ I 



- Tr In — - i Tr In V — Tr ( A-^TB-^V 

2 Bo ^i2n V 



Ao := 



:p=p=0; 



j=p=p=0 



(2.10) 



Recalling that A^^, are the scalar and fermion matrix propagators in the presence of 
external sources, the one-loop functional Wl=i is seen to be a sum of the bosonic one-loop 
functional | Tr \n{A/ Aq), the fermion-loop functional —i Tr hi{B / Bq) and a mixed one- 
loop functional where scalar and fermion propagators alternate. In order to determine 
the ultraviolet divergences that occur in the mixed term in ( |2.1(j| ), the calculational 
inconveniences mentioned in Sect. |l| are encountered. 

These problems can be circumvented by reorganizing the three parts of Wl=i into a 
more compact form, using the notion of supermatrices, supertraces, etc. (see for instance 
P, P]). Combining the bosonic and fermionic fluctuation variables in a multicomponent 
field 



A 



V 

—T 

rj 



(2.11) 



m 



can be written as 



(2.12) 



The explicit form of the supermatrix operator K follows immediately from the second 
line in (^): 

A r -r^ 
-r"^ -5^ 
r 



K 



(2.13) 



The one-loop functional of connected Green functions can now be written in compact 
form P in terms of a supertrace 



L=l 



(2.14) 
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With the notation 



Str O = J d'^x str(x|0|x) 



I distinguish supertraces with and without space-time integration. 

For actual calculations, the form of the supermatrix operator K defined in ( p.l3|) is 
not the most convenient one. Applying a similarity transformation to K, the generating 
functional can also be written as 



with 



K' 



L=l 



A 



(2.15) 



r " 

fxB . (2.16) 

fiB^ 

The arbitrary mass parameter /i introduced in ( p.l6| ) guarantees equal dimensions for 
all entries in K' {[K'] = [A] = 2). Although this quantity does not, of course, appear 
in any final result, it turns out to be quite helpful for the inspection of expressions at 
intermediate stages of calculations. 

In the proper-time formulation, the one-loop functional assumes the form 



Wl=i 



r^Str (e 
2 Jo T 

dr 



itK' _ ^irK'o 



2 JO T J 



(2.17) 



which is just the desired super-heat-kernel representation. Note that the convergence of 
the integral at the upper end (r —>■ 00) is guaranteed by the small imaginary parts present 
in the bosonic and fermionic differential operators A and B, which are ensuring at the 
same time the usual Feynman boundary conditions. (For a free theory A = — □ — -|-ie, 
B = i P — m + ie.) On the other hand, the behaviour of the integral at the lower end 
exhibits the divergence structure of the theory under investigation. 

As long as we are only interested in those parts of the one-loop functional that are at 
most bilinear in fermion fields, the supermatrix K' can be reduced to the simpler form 



K" 



A F 
2J1T fiB 



(2.18) 



such that the one-loop functional reads 



- Str In — - — - Tr In h . . . 

2 K'^ 2 Bo 



(2.19) 



The terms omitted are at least quartic in the fermion fields. 
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3 Scalars Interacting with Heavy Fermions 



In the case of chiral perturbation theory with heavy baryons, the fluctuation action ( |2.9|) 
generated by the lowest order meson-baryon Lagrangian in the mesonic and 0{p) 

in the baryonic part) has the general form 

^(2) = -^fiD.D'^ + Y)^ + fj{a + (3,D^)i + ^^(5 - ^^P^)r/ + l)tv,V^^ii , (3.1) 
where 

J = a-V;p^, (3.2) 
= 1 , y . = . 

X^, Y and are bosonic (matrix-) fields, whereas a and are fermionic objects. The 
form of 5 in ( p.2| ) is required by the reality of ( |3.1| ). is the usual velocity vector 
introduced in the heavy mass expansion. Apart from the condition w ■ /3 = (which 
is indeed satisfied in heavy baryon chiral perturbation theory), no further assumption 
about the terms entering in ( p.l| ) is made in this section. The discussion will therefore 
apply to a rather large class of theories of scalars interacting with heavy fermions, not 
necessarily related to chiral perturbation theory. The further specialization to the pion- 
nucleon system is reserved until the next chapter. 

The action (p.l|) is invariant under local gauge transformations 

i{x) R{x)i{x) , R{xfR{x) = 1 , 

rj{x) U{x)rj{x) , U{xyU{x) = 1 , 

Xfj, RdnR~^ + RX^R^^ , 

Y RYR~^ , (3.3) 

U ^ ud^u-' + uf,u-^ , 

a UaR^^ , 

^ UP^R-' . 

Consequently, also the divergent part of the one-loop functional exhibits this symmetry 
property UTof. The matrix- fields Y, a, (3^ together with their covariant derivatives 



<9/,a + f^a - aX^, 
d^(3u + - (3.4) 



and the associated "field-strength" tensors 

Xfj_i, := d^Xy — dyXfj_ + [X^, X^], 

fp^y ■■= d^fy - dyf^ + [/^, fy] (3.5) 

are therefore the appropriate building blocks for the construction of a gauge-invariant 
action. 
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The general heat-kernel formalism of the preceding section will now be applied to 
In this case, the matrix-operators A, B, T and F defined in ( p.9|) are given by 



A = -D^-Y, B = tvV, T = a + (3-D, T = 6 - (3 ■ V . (3.6) 
As I am considering only terms at most bilinear in the fermionic variables, the form ( |2.18| ) 



for the supermatrix operator is the appropriate one. Employing the method of Ball |jTT 
the relevant diagonal space-time matrix element can be written as 



str (x\e 



rK" 



\x) 



str 



str 



d k{x\e 



rK'' 



\k) {k\x) = str 



d'^k 
(27r)« 



ikx itK" —ikx 

e e e 



(3.7) 



with 



K" 



(3.8) 



-D2 _ y + fc2 + 2?fc ■ D -p-V + ik-P) 

v^(a + (3- D-ik- (3) n{iv -V + vk) 

The further evaluation of this expression is considerably simplified by the observation 



that in the following intermediate steps we may restrict ourselves to constant fields |Tl 
X^, a, Pf^, f^, Y = —X^. As the final result for the one- loop divergences has to be gauge- 
invariant, no information is lost and the full expression for space-time dependent fields 
is recovered by the substitutions 



-X^ - 


- Y, 


-[x„x'] - 








f^a - aX^ - 




ffiPu — PuX^ - 




\f fv] ~ 





In this approach, (|3.71 ) reduces to the much simpler expression 

d'^k 



str (x\e 



irK" 



\x) 



str e^^+^ 



(3.9) 



(3.10) 



with 



M 



N 



IT 



k'^ + 2ik-X 

ii{iv ■ f + V ■ k) 

6 -P - f + ik-p 

a + p-X -ik- P 



(3.11) 



Let us first consider the part bilinear in the fermionic matrix (generating the terms 
of the form a . . .a,a . . . P^. . .a, P^. . . P^). The corresponding part of the generating 
functional (|2.1CI|) is just 

■^^ilp ^ := -I Tr {A-^VB-^V) . (3.12) 
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The appropriate decomposition of the exponential in (|3.1CI|) can be performed by using 
Feynman's "disentanghng" theorem [[12| : 



exp(M + A^) = exp M P, exp / ds N{s) 

Jo 



(3.13) 



with 



N{s) := e-'^'Ne'^ 



and 



n=0 



Fs exp dsN{s):=J2 dsi ds2 . . . dsn N{si)N{s2) . . . N{sn) ■ (3.14) 



(In the mathematical literature, (|3.13| ) is also known as "Duhamel's formula".) Picking 
out the part bilinear in A^, 



str e*^+^ = ds ds' str U^-s)m j^^is-s')M j^^s'm' 
Jo Jo ^ 



+ ... 



(3.15) 



a few simple manipulations lead to 



str e 



M+N 



-2iiT^ / dz e*-[(i-^)'='+^^''-^l tr \{5-p- f + ip- fc)e-""^^-^ 

Jo L 

{a + f3-X-if3- k)e-^^^^-'^''-^] + ... 



(3.16) 



After integration over z, the /i-dependent terms cancel once the proper-time and the 
momentum-space integrals are applied. The remaining contribution to Wl=i assumes the 
form 



WL=i\r 



r...r 



d<ix r-e-'' I T^e»-' tr 



t 



(2 



(6 - p- f + il3-l/t)e 



-tvf 



{a + p-X-ip- l/t){f + 2itl ■ Xy^ 



(3.17) 



where a suitable change of the integration variables has been performed. The divergent 
part (for d — > 4) can now be easily isolated: 



W^" I- 



1^ 



tr \ [5-p.f)vf{a + P-X) + ^^P-l{vffP-l 



+ 



{8-P-f){a + P-X) + i{5-P-f)vfP-l 



-iP-lvf{a + P-X) + j^P-l{vffp-l 



2il-X 



+ 

-p-lp-l 



i{6 - p ■ f) p ■ I + ip ■ I {a + p ■ X) - p ■ Iv f p ■ I 
«z(/-X)n 



2\3 



A{l-Xy 
(3.18) 



The necessary formulas for the /-integration are given in the Appendix. In the last step, 
one has to identify the appropriate gauge-invariant combinations (constituting a non- 
trivial check of the calculation) and reconstruct the full result by using ( p.9|) . In this way, 
I finally obtain: 



W 



div |_ 

L=ilr...r 



487r2(rf-4) 



d^x tr [-I2av ■ Va + 6 \a(3^X^"'v^ + [i^aX^'^v^ 



+6 ^^{v ■ V/5,)X^'^ + 4 V ■ VX^" + 2 (3^(3,V^'X'"vp] . 



(3.19) 



Note that ( p.l9 ) has to be real, which is another independent check of the result. 

The remaining part of the generating functional with the fermionic operators F, F 
turned off, 



L=ilr=r=o 



— Tr In 2 Tr In — 

2 Ao Bo 



(3.20) 



does not require any additional effort. A simple calculation (involving a Gaussian 
momentum-space integration) gives 

1 



Tr InALi 



div 



(47r)2(d-4) 



d^x tr f ^X^^X^-" + ^Y' 



(3.21) 



which is the standard result obtained by 't Hooft [|T^ using diagrammatic methods. 

The second term in (|3.2CI| ) vanishes identically, as it corresponds to the closed loop of 
a "light" fermion component in the heavy mass expansion: 



Tr In 5 



d'^x 



d'^x 



dt 
T 
°° dt 
T 



-t' 



d'^k 

dH 



tr /'g«*(™-^+^-'=)]_ 



(2vr) 



which follows from 



dH 



Av-l 



. 



4 Renormalization of the Pion— Nucleon Interaction 



The functionals ( p.l9| ) and ( p. 21] ) are the basic formulas for the analysis of the one-loop 
divergences to 0{p^) in heavy baryon chiral perturbation theory. They can be applied to 
both the two-flavour and the three-flavour case. In the following I shall confine myself to 
chiral SU{2). 

The starting point for the formulation of the effective field theory of the pion-nucleon 



system is QCD with the two light flavours u, d coupled to external Hermitian fields |13 



r - 



u 
d 



(4.1) 
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Cqqy) is the QCD Lagrangian with m„ = rud = 0, S and P are general two-dimensional 
matrix fields, the isotriplet vector and axial-vector fields Vfj_,A^ are traceless and the 
isosinglet vector field is included to generate the electromagnetic current. 

Explicit chiral symmetry breaking is built in by setting S = quark = diag [mu,md]. 
The chiral group G = SU{2)i x SU{2)r is spontaneously broken to the isospin group 
SU{2)v- It is realized non-linearly on the Goldstone pion fields 0: 



ul{4>) 



^ gLUL{(p)h{g,(f)) \ 
gRUni(l))h{g,(f))-\ 



g = {gL^gn) e G 



(4.2) 



where ul.ur are elements of the chiral coset space SU{2)l x SU{2)r/ SU{2)v and the 
compensator field h{g, 0) is in SU{2)v. 
The nucleon doublet \E' transforms as 



P 

n 



^' = h{g,cl>)^> 



(4.3) 



under chiral transformations. The local nature of this transformation requires a connec- 
tion 



ul 



V„ - A,. 



(4.4) 
(4.5) 



^ 2 

in the presence of external gauge fields 
to define a covariant derivative 

V^vl> = (9^ + r^-,V;)vl>. (4.6) 
To lowest order in the chiral expansion the effective Lagrangian of the pion-nucleon 



system is [ffl, 13 



with 



-eff 



u. 



gA 



(4.7) 



X 



2B{S + iP), 



X+ 



F, m, gA are the pion decay constant, the nucleon mass and the neutron decay constant in 
the chiral limit, whereas B is related to the quark condensate. (. . .) stands for the trace 
in flavour space. 

The heavy baryon mass expansion of (|4.7| ) is obtained by introducing velocity- 
dependent flelds 



(4. 
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leading to 



F2 
4 



jO-es = —{u^u^' + x+) + Ny{zvV + gAS-u)N^ + ... 



(4.9) 



The additional terms involving the "heavy" fermion components are irrelevant for our 
present purposes. (For a more detailed discussion the reader is referred to |15l-) The 
only dependence on Dirac matrices in (^4.9|) is through the spin-vector matrices 



S-v = 0, S' = --l 



(4.10) 



which obey the (anti-) commutation relations 



{S^^ S"} = -[v^v"" - g''") , [S", S"] = ie^''P"vpS, . 



(4.11) 



To obtain the associated second-order fluctuation Lagrangian ( [4.9|) is expanded 
around the classical fields 0ci, N^^ci- In the standard "gauge" Ui?(0ci) = ^l(0ci) =: u{(j)ci) 
a convenient choice of the bosonic fluctuation variables (i = 1, 2, 3) is given by |l! 



ur{<P) = M(0d) e^^, = u\(l)d) e e (0ci) = , (4.12) 



where r denotes the Pauli matrices. For the fermion fields 1 write Ny = N^^ci + V- In this 
way I get 



(2) 
eflf 



-[v-u, Ti]^i - gASf,Ti{d^jQ 



1_ 

+rj{iv ■ V + gAS ■ u)r] 



V 



-[v-u, - gASf,Ti{dfjQ) 



(4.13) 



where 



(4.14) 



Note that the quantities A^^, u^, F^, x+ in ( [4.13| ) are to be taken at the solutions of the 
classical equations of motion. (The subscript "cl" has only been dropped for simplicity.) 
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It is now easy to verify that the action associated with ( [4.13|) can indeed be written 
in the standard form (|3.1j ) by setting 





= + 9tM , 9ij = 


Y 


= a + s , Sij = 




= - iV^ - iVf^QAS 











8F2 

9A 



N,[Ti,Tj]N, 



1,2,3, 



4F2 



A^^ {2 6ij S ■ u — Ti S ■ u Tj — Tj S ■ u Ti) 



a = 1,2 , 



(4.15) 



Let us first consider the one-loop divergences generated by (|3^) . Using (|4T5|) , 



tr 



tr cr^ + 2 tr (as) + • • • 



and 



where 



tr {X,,Xn = tr 7^.7^^^ + 4 tr {^..{d'^g'' + [7^ g'])} + • • • 



(4.16) 



(4.17) 



7^1' := df,-fu - d^-f^ + [7^, 7^] . (4.18) 

The first terms on the right-hand sides of (|4.16|) and (|4.17|) are purely mesonic; they 
determine the divergence structure of the well-known Gasser-Leutwyler functional of 
0{p^) [|l^]. The second ones are bilinear in the fermion fields, whereas the dots refer to 
irrelevant terms ~ {N^ . . . Ny)"^. To facilitate the comparison with I write the fermion 
bilinears extracted from ( 3.21| ) in the following form: 



W] 



div I _ 

L=ilr=r=o 



div 



1 



87r2F2(d-4) 



El . 



(4.19) 



For Si I find 



Si 



where 



(4.20) 
(4.21) 

This result agrees with the corresponding expression in (36) of Note that I have used 
several SU{2) relations to arrive at a simpler form for Si in (|4.20| ). 

The one-loop divergences originating from ( |3.19D are again presented in the form 



div I _ 

L=ilr...r 



d^x N,, S^'" Ny 



-idiv 



87r2F2(d-4)^ 



^2 • 



(4.22) 



11 



Inserting ( |4.15| ) in ( p.l9|) , I obtain: 



i |-[2(i; ■ uY + {{v ■ u)'^)]v ■ V + -f ■ u{v ■ Vv ■ u) + -(f ■ u{v ■ Vf ■ u)) 
+ Qa ■ ^(-S* -uv -u) + ^(S* ■ u{v ■ uf) - S^v''\rf,^,v ■ u] 



+ w\ ■ V)' - - {V^V.^v") + te^'''''vpS^[2T^,v ■ V + (i; ■ Vr^.)] 

_ A . di^Au ■ u + 3x+)) - ^(4m ■ M + ?,x+)v ■ v| 

+ g\ ■ ■ V)=^ - 2^^(^ ■uV^;)S'' - ]^{v ■VS-u)vy 

((t; ■ VfS ■ u) - ^u,{u^S ■ - ^5 ■ 
+ z^^S"^ |[2(5 • uf - 4{{S ■ uf)]v ■V+'^{v-VS-u)S-u 
+ ■ u{v -VS-u)- A{S ■ u{v ■ VS ■ 

+ 9'AS,{l{S-uf-^{{S-uf)]s'^ , (4.23) 

which is in agreement with the corresponding result in (Note that "H — " in the fourth 
hne of (53) in p should be read as a minus sign.) 

5 Conclusions 

I have shown that the super-heat-kernel technique constitutes the appropriate theoretical 
tool for analyzing the one-loop divergences in systems with (non-renormalizable) boson- 
fermion interactions. I recall here the essential ingredients that were combined to arrive 
at an efficient computational scheme: 

• The one-loop functional is written in terms of the superdeterminant of a suitably 
chosen supermatrix operator. 

• The associated super-heat-kernel representation is the appropriate form of the one- 
loop functional for studying its divergence structure. 

• It is easier to determine the diagonal heat-kernel matrix elements directly by in- 
serting a complete set of plane waves instead of calculating the Seeley-DeWitt co- 
efficients with two different space-time arguments and taking the coincidence-limit 
at the end. 

• The heat-kernel-representation is perfectly well defined also for supermatrices with 
first-order (fermion) differential operators^]. 



^Note, however, that "squaring" of the fermionic differential operator may simplify the analysis in 
theories where the full relativistic Dirac operator is still present 0|. 
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• The second-order fluctuation action is invariant under a local gauge transforma- 
tion. As a consequence, this symmetry property is also shared by the divergence 
functional. 

• At intermediate stages, the calculation can be carried out with constant (classical) 
fields, avoiding cumbersome manipulations with derivatives acting on space-time 
dependent objects. At the end, the general result is recovered by gauge invariance. 

• Feynman's disentangling theorem allows the proper decomposition of the exponen- 
tial of a sum of non-commuting terms. 

• With the divergence functional given in compact form, the one-loop renormalization 
of effective quantum field theories becomes an easy task, requiring only a few purely 
algebraic operations. 

The application to heavy baryon chiral perturbation theory with two flavours served as 
an explicit example. A previous result for the counterterms to 0{p^) was confirmed. 

With the super-heat-kernel method at hand, the systematic study of effective field 
theories at the one-loop level is simplified considerably. I am giving here a small selection 
of possible applications and extensions of the present work: 

• The treatment of the meson-baryon interaction with three flavours is completely 
analogous to the two-flavour case disussed before. 

• The inclusion of fields with higher spin (photon, A-resonance, etc.) is straightfor- 
ward. Their components are simply added to the bosonic and fermionic sectors, 
respectively. 

• The completion of the one-loop renormalization for the pion-nucleon interaction 
up to 0{p'^) may be achieved by a suitable extension of (|3.1|). 



For the analysis of fermionic bound states, the complete form ( |2.16|) of the super 



matrix operator must be used, as terms quartic in the fermion fields are relevant in 
this case. 



In analogy to the mesonic case [^, the super-heat-kernel representation might also 
be useful for the finite part of the one-loop functional with two external baryons. 



Acknowledgements 

I am indebted to Gerhard Ecker, Jiirg Gasser, Joachim Kambor and Marc Knecht for 
helpful discussions and useful comments. 



13 



Appendix 

I consider first the integrals 



with an arbitrary four- vector f^. The fn{d) are given by 



w-i r(rf-2) 



fJd) = i-iY-' )r-i — (A.2) 



and 



The momentum space integrals occurring in ( |3.18|) are now obtained by differentiating 
( [A.l| ) a sufficient number of times with respect to and setting f ^ = 1 at the end: 



a L e ^ 



{2'kY P 47r2 



(A.4) 



If^L — > 7r^{9t^u-^v^v^) , (A.5) 



{27rY P 27r2 
d^^l e"'-' 



^ (A.6) 



{27tY {Pf ^ 87r2 ' 



dwM- ^ A (^7/.. - 2t;;.t^^) , (A.7) 



{2ttY {Py ^^^"^^ ^ J^[-(9,.uVp + ...)+Av^v^Vp] , (A.8) 

e*''^ 1 

^ ^[-(^7M-^^P + ...) + 2t;;.t;.t;p] , (A.9) 



+8vpV,yVpVa] , (A. 10) 



d'^l e™-' , , , , , 1 



IplvlplJr > 2 [~ {dtiudpaVr + ...)+ 2 {QpvVpV^Vr + • • •) 



{2'kY {Py ' " 1927r2 

-8t;^t;i,t;pt;^t;^] . (A. 11) 

The dots indicate the necessary symmetrizations. 
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